Spiky density of states in large complex Al-Mn phases 
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First-principle electronic structure calculations have been performed in crystalline complex phases 
^-Al4.i2Mn and A-AU.eMn using the TB-LMTO method. These atomic structures, related to qua- 
sicrystalline structures, contain about 560 atoms in a large hexagonal unit cell. One of the main 
characteristic of their density of states is the presence of fine peaks the so-called "spiky structure" . 
From multiple-scattering calculations in real space, we show that these fine peaks are not artifacts 
in ab-initio calculations, since they result from a specific localization of electrons by atomic clusters 
of different length scales. 



I. INTRODUCTION 



Al based Quasicrystals, approximants and related com- 
plex metallic phases reveal very unusual transport prop- 
erties (for review see Ref. [l| and Refs. therein). Indeed, 
they exhibit a very high resistivity although the den- 
sity of states (DOS) at the Fermi energy Ep is not very 
small. This means that their transport properties are 
due mainly to a small difFusivity of electrons. Recently 
, we proposed a new mechanism for the conductivity in 
those systems where the velocity of charge carriers is very 
low. This mechanism is related with a specific localiza- 
tion of electrons by complex atomic structures. Several 
numerical studies argue that this localization leads to 
fine peaks in the DOS, called "spiky peaks". Thus, spiky 
DOS were predicted in icosahedral small approximants 
(for instance aAl-Mn-Si l/l-approximant [3, 1/1 Al- 
Cu-Fe (3, 1/1 Al-Pd-Mn [|). In the case of crystalline 
approximants, it is obvious that these peaks are associ- 
ated to a small electron velocity (flat dispersion relations) 
Q. Such properties are not specific of quasicrystals as 
they are also observed in many crystal related to qua- 
sicrystals, therefore it does not come from the long range 
quasiperiodicity. Indeed, we have predicted that fine 
peaks in the DOS could come from electron confinement 
by atomic clusters characteristic of the quasiperiodicity. 
This is not in contradiction with a Hume-Rothery mech- 
anism for stabilizing these phases because this tendency 
to localization has a small effect on the band energy 
Q. The existence of spiky DOS in quasicrystals and 
approximants with large unit cell is however still much 
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debated experimentally (Refs. 

therein) and theoretically (Refs ^ ^ 

therein). Very recently, R. Widmer et al. [l3|, found a 
clear signature of spiky local DOS near Fermi level in 
scanning tunneling microscopy and spectroscopy at low 
temperature. 

In this paper, we present ab-initio calculations and 
multiple scattering calculations that prove the existence 
on fine peaks in the DOS of complex Al-Mn crystals with 
large unit cell. 




II. FIRST-PRINCIPLE STUDY 

The /i-Al4.i2Mn phase is crystallized in a hexagonal 
structure (PGa/mmc space group) with large unit cell: 
a = 19.98(1)1 and c = 24.673(4) A [H. Each unit cell 
contains 110 Mn atoms and ~453 Al atoms. As explained 
in Ref. [l3|, we have modified slightly the Wyckoff po- 
sition of one Al atom in order to avoid too short in- 
teratomic distances. The A-AU.gMn crystalizes also in 
PGa/mmc with a = 23.382(9) A and c = 12.389(2) A 
[ISj . Its atomic arrangements have many similarities with 
those of /i-Al4.i2Mn. 

DOSs are calculated by using the self-consistent 
tight-binding Hnear muffin tin orbital (TB-LMTO) 
method p^ . in the atomic sphere approximation (ASA) 
without spin polarization. The LMTO basis includes all 
angular momenta up to Z = 2. The Brillouin zone (BZ) 
integration on the k space was done by the tetrahedron 
method with Nk A:-points in a BZ. LMTO DOS of /i- 
Al4.i2Mn have been already calculated by D. Nguyen- 
Manh and G. Trambly de Laissardiere [l3] for analyzing 
the origin of rnagnetism in complex aluminides related to 
quasicrystals [l3,[3- Here, we focus on the shape of the 
DOS without spin-polarization. To increase the accuracy 
of calculation, we increase the number N^- In Ref. 
TB-LMTO calculation in /i phase was performed with 
Nk = 9, whereas Nk = 3375 in actual calculation for fi 
phase, and Nk = 1458, for A phase. 

The total DOSs in /i and A phases are very similar 
(figures [T] and [2]) . At low energies, the parabola due to 
the Al nearly-free states is clearly seen. The large d band 
(figure Hfb)) from —2 up to leV is due to a strong sp- 
d hybridization [1, d, [11, [lO]. The sum of local DOSs 
on Al atoms, shown in figure (H^c), is mainly sp DOS. 
As expected for a Hume-Rothery stabilization, there is a 
wide pseudogap in sp DOS near Ep 0, H, H, [ll, [H, HO] . 
The pseudogap in the total DOS is narrower because the 
d states of Mn atoms must fill up it partially fl^. As 
shown on the insert figure [U DOS around Ep consists in 
a set of spiky peaks. Nevertheless the resolution is not 
well define in the TB-LMTO calculations [4], therefore 
these calculations are not enough to prove the existence 
of spiky peaks. 
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FIG. 1: Total LMTO DOS in /i-Al4.i2Mn. Insert: Total DOS 
around E — 0. 
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FIG. 2: LMTO DOS in A-AU.eMn: (a) Total DOS, (b) local 
Mn DOS and (c) local Al DOS. 



III. SPIKY DENSITY OF STATES ? 

To prove the existence or not of fine peaks in DOS, we 
calculate the DOS for a large assembly of atoms by using 
a multiple scattering approach in real space. Therefore 
there is no error due to Brillouin zone (BZ) integration as 
with TB-LMTO method. Multiple scattering calculation 
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FIG. 3: Variation An(£;) of the total DOS due to X Mn 
atoms in free electron matrix. Mn atoms are located on the 
Mn sub-lattice of /i-Al4.i2Mn. (a) AT = 1 is the Virtual Bound 
States, (b) X = 110 corresponds to 1 unit cell of ^ phase, (c) 
X = 880, 8 unit cells of ^ phase, (d) X = 1980, 18 unit cells 
of /^ phase. 



is not an ab-initio calculation, but it allows to discuss 
qualitatively the existence of spiky peaks. 

In a standard description of intermetallic alloys, one 
starts from a muffin-tin potential, which varies in spheres 
centered on the atoms and is constant in the space be- 
tween the spheres. As far as the band structure is con- 
cerned an atom is entirely characterized by its scattering 
properties for incident plane waves. The effect of a group 
of atoms on band structure will be characterized by its 
T-matrix, that is by its scattering properties of plane 
waves (Ref. Q and Refs. therein). To simulate Al(Si)- 
Mn phases, we consider only the potential due to the Mn 
atoms, which are strong scatterers, and we neglect the 
potential of all the Al (and Si) atoms that are consid- 
ered as weak scatterers. Furthermore, we consider that 
the Mn d orbitals are the only source of scattering for the 
conduction electrons (i.e. the sp electrons) . Within these 
approximations a complex Al(Si)-Mn phase becomes an 
assembly of X atoms of Mn in a jellium simulating Al 
(Si) atoms (free sp electron matrix). The total DOS of 
this system is: 



(1) 



where is the DOS of free electrons without sp-d 
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FIG. 4: Variation An{E) of the total DOS due to X Mn atoms 
in free electron matrix. Mn atoms are located on the Mn sub- 
lattice of a-Al-Mn-Si cubic 1/1 approximant. (a) X = 1 is 
the Virtual Bound States, (b) X = 12, 1 Mn icosahedron. 
(c) X = 192, 16 Mn icosahedra on bcc lattice, (d) X = 648, 
54 Mn icosahedra on bcc lattice, (e) X = 1536, 128 Mn 
icosahedra on bcc lattice. 



interaction, and An, the variation of DOS due to the Mn 
atoms. An{E) is calculated by using the Lloyd formula 
0) im , which is a generalization for several Mn atoms of 
the Friedel sum rule for one Mn atom in jellium: 

An{E) = -^Tr (miE)-' ^MiE)^ . (2) 

The matrix elements of M(i?) are calculated form the 
T-matrix: 



M 



G 



(3) 



where G is the propagator that coupled d orbitals of 
different Mn scatters, and 



(4) 



tx is the element of the transfer matrix due to the Mn 
scatter indexed by x. depends only on the phase shift 
Sd in wavefunction of the sp states due to their scattering: 



tan Sd{E) 



with 



E-E, 



2T^TTV^n{;^'''(EF). 



(5) 
(6) 



where 2r is the width of the d resonance and Ed the 
on-site d energy. V is the coupling between one d orbital 
and sp eigenstates [l^. When X = 1 (impurity limit), 
the total sp DOS, Ugp, is not modified by the sp-d hy- 
bridization, and the variation An of the total DOS due 
to the Mn atom is a Lorentzian (Virtual Bound State) 
centered on energy Ed and with a width equals to 2r. 

To simulate complex Al(Si)-Mn, the only parameters 
are the values of Ed and T, and positions of Mn atoms 
in the real space. We use the following realistic values 
Ed = 10.9 eV and 2r = 3 eV, and atomic po- 
sitions determined from X-ray refinements. The DOS 
of several assemblies of Mn atoms are thus determined 
for Mn belonging to the Mn-sublattice of /^-phase (fig- 
ure [3]) and a-AlMnSi 1/1-approximant (figured]). For fi 
and a structures, fine peaks appear in the DOS when 
a large number of Mn scatterers is taken into account 
(large value oi X). Widths of these peaks are about 
40 — 150 meV which agrees with ab-initio calculations 
(figures [T] and [2] and Ref. [11] )■ Each fine peak is a 
signature of states confined by atomic clusters, so-called 
"Cluster Virtual Bound states" When the number 
of Mn scatters increases, new peaks appear. This means 
that atomic structure confines electrons on several length 
scales corresponding to different sizes of clusters, respec- 
tively. 

Summarizing, multiple-scattering calculations have 
been performed in real space to simulate DOS of com- 
plex Al(rich)-Mn phases related to quasicrystals. They 
show that a localization of electrons by large assemblies 
of atoms leading to fine peaks in the DOS. This confirms 
the existence of the spiky DOS found by ab-initio calcu- 
lations. 

This work owes much to discusion with D. Mayou and 
D. Nguyen-Manh with whom ab-initio study of /x-phase 
has been initiated. The computations have been per- 
formed at the Service Informatique Recherche (S.I.R.), 
Universite de Cergy-Pontoise. I thank also Y. Costes 
and B. Marir, S.I.R., for computing assistance. 
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